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Abstract 

Quantum phase transition in dimerized antiferromagnetic Heisenberg spin chain has 
been studied. A staircase structure in the variation of concurrence within strongly coupled 
pairs with that of external magnetic field has been observed indicating multiple critical 
points. Emergence of entanglement due to external magnetic field or magnetic entangle¬ 
ment is observed for weakly coupled spin pairs in the same dimer chain. Though closed 
dimerized isotropic XXX Heisenberg chains with different dimer strengths were mainly 
explored, analogous studies on open chains as well as closed anisotropic (XX interaction) 
chains with tilted external magnetic field have also been studied. 


1 Introduction 


Quantum phase transition (QPT) is a phenomenon where the non-analyticity in the energy level 
gives rise to a situation where energy fluctuation at the critical point vanishes [l]. Entanglement 
in a quantum system may be considered as a witness of this phenomenon (see e.g. |2,3|). On the 
other hand, entanglement being an important resource 00 in quantum information processing, 
this phenomenon is important from the perspective of its physical realization too. 

Spin-1/2 chains are some of the simplest qubit systems considered in the literature (see 
e.g. [sqTl] ) for the implementation quantum processor. Again the case of a two spin system, 
interacting through an isotropic XXX Heisenberg Hamiltonian in presence of an external mag¬ 


netic field, the entanglement of formation has been calculated by Nielsen 12 , which was further 


extended for a spin chain under periodic boundary conditions with a systematic investigations 
on the variation of entanglement between various sites due to that of an external magnetic field 
and temperature (see HD- The pairwise entanglement vanishes even at a vanishingly small 
temperature if the magnetic field exceeds a critical value which characterises quantum phase 
transition. The effect of magnetic entanglement on isotropic XY model and thermal entan¬ 


glement for anisotropic XY model 13 and with non zero external magnetic field the effect of 


anisotropy on the same system had also been studied 14 . The dynamical phase transition in 


anisotropic XY model was also been studied 15 


The effect of non uniform magnetic field on thermal entanglement for the two qubit system 
with isotropic Heisenberg interaction [16] and for a two-qubit Ising model [l7| were studied. 
Gong and Su 1181 have investigated pairwise thermal entanglement in S = \XY chain (including 
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its thermodynamic limit N —» oo) and compared their results, with mean field frame work. 
They have found that in the case of first model there exist a common T c , independent of B for 
both ferromagnetic and anti-ferromagnetic interactions. It had also been found that in such a 
spin system Tq is independent of B for AF interactions. The studies were further extended for 
Heisenberg model with next to nearest neighbour interaction 
of QPT was noticed. Other possibilities like existence of tilted external field 
with spin-1/2 and spin-1 elements [21], were also explored. Far from equilibrium QPT has been 


and existence of two types 
mixed chain 


20 


observed in an open Heisenberg XY spin-1 chain characterized by a sudden appearance of long- 
range magnetic order in the non-equilibrium steady state as the magnetic field is reduced 22 


In fact for last one and half decade or so, QPT in Ising chains and Heisenberg chains was 
extensively studied by different workers, only some of which we mentioned here. 

Though the dimerized versions of Heisenberg spin chains were studied since a long time back, 
they got renewed interests from the perspective of quantum communication. On the other hand, 
for quite a long time spin chains are proposed candidates for communication medium between 

and references therein). But very recently a scheme of 
using dimerized spin-1/2 XXX model. Another 


two quantum processors (see e.g. 23 


information transfer has been proposed 24 


scheme to generate end to end entanglement has also been proposed 125] using dimerized XXX 
and XXZ model. Although phase transition in dimerized spin-1/2 chain was studied in connec¬ 
tion with transition betweeen gapped and ungapped energy spectra (see e.g. 26,27]), quantum 


phase transition for the dimerized cases using thermal and magnetic entanglement has not 
been studied yet according to the authors’ knowledge. The present paper explores this field for 
spin-1/2 dimerized Heisenberg chain. 


2 Formalism 

The exchange Hamiltonian of a Heisenberg spin chain may be generalised as 

He = X^X+i + Vf^+i + (1) 

i 

where J x , J y and J z are the coupling strengths along three orthogonal directions. In addition if 
some external magnetic field B is applied on the system, there will be additional on-site energy 
terms, generally called as Zeeman energy terms, as 

H B = Y J B.a i . ( 2 ) 
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Figure 1: Schematic diagram of a dimer chain 

The dimerized Hamiltonian for a chain of N spins in presence of an external uniform mag¬ 
netic field acting along z-axis is given by 


N 


N 


h = +j£> + (-y+U]*. 




(3) 
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are the Pauli spin matrices and the dimer strength is denoted by <5; with 0 < 5 < 1. It is quite 
obvious from the expression that the coupling strength becomes strong and weak alternately 
for consecutive pairs which characterises a dimer. For the AFM chains the interaction strength 
J > 0 and for the FM cases J < 0. For a closed chain i runs from 1 to N in the second 
summation within the Hamiltonian expression, eq. (J3|, with the boundary condition epv+i = o \; 
while for open chain i simply runs from 1 to IV — 1. For an external magnetic held tilted by an 
angle 6 with £-axis, eq. ([3]) may be modified as 

N N N 

H = B cos 9 of + B sin 6 erf + J ^^[1 + (—1 ) l+1 h]a i .(T i+ i. (4) 

2=1 2=1 2=1 

The system is taken to be in thermal state i.e., the state of the system when it is at thermal 
equilibrium with a thermal bath (environment) at temperature T, such that the density matrix, 
p{T ) = exp (—^)/Z where Z is the partition function and k is Boltzman’s constant. In this 
communication we are reporting the variation of pairwise thermal entanglement for the strongly 
coupled and weakly coupled spin pairs separately with that of magnetic held as well as that of 
dimer strength 6 for a closed AFM chain(along with the open ended ones) of length N. 

For this purpose hrst we need to find out the pairwise entanglement of qubit pairs. This 
may be obtained from the reduced density matrix for that pair of qubit, evaluated by tracing 
out the contribution of other qubits (or in other words taking partial trace on those degrees of 
freedom) from the total density matrix p{T) of the system. Thus the reduced density matrix 
of the subsystem (spin pair) of our concern may be expressed as 

Pr = Tib \p(T )\, 


where the complementary part of the system is named as subsystem B. 

Now in general p^ will be of mixed state. For such a system the suitable measure of 
entanglement may be considered as entanglement of formation which may be expressed, for a 
two qubit system, as a monotonic function of a simpler variable named as concurrence 
From reduced density matrix the following product matrix was calculated 


R = PR^y ® °y)P*R ( K (7 y ® Vy) ( 5 ) 

If we represent Ads with i = 1,2, 3,4 as the square root of the eigenvalues of R in descending 
order then concurrence can be expressed as 


C = max(0, Ai — A2 — A3 — A4) (6) 

In this method we have described the standard bases,100), |01), 110), 111) must be used. In 
this study we have calculated the concurrence indicating the thermal entanglement between 
different spin pairs. 


3 Multiple Transitions and Staircase Structure 

Following the same line of investigation as done in the ref. |2|, we invetigated the pairwise 
entanglement of nearest neighbour pairs for a closed chain as quantified by the concurrence; 
but in this case for two types of pairs viz., strongly coupled and weakly coupled separately. 

It is quite apparent, from the £gure^2] and [3] , that in the case of a dimer too, there is a sharp 
critical value of B = B c across which the pairwise entanglement, both in the case of strongly 
coupled as well as of weakly coupled pairs abruptly vanishes. At very low temperatures the it 
is more abrupt due to the minimal presence of thermal fluctuations indicating some quantum 
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Figure 2: Variation of entanglement between strongly coupled spin pairs with temperature 
and external magnetic field for closed XXX chain. N = 12, 5 = 0.2, 9 = 0 



Figure 3: Variation of entanglement between weakly coupled spin pairs with temperature and 
external magnetic field for closed XXX chain. N = 12, 5 = 0.2, 9 = 0. 



Figure 4: Variation of entanglement within strongly coupled pair with magnetic field at 
constant temperature i.e with fixed ksT = 0.1. N = 12,6 = 0.2, 9 = 0. 
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Figure 5: Variation of entanglement within strongly coupled pair with magnetic field at 
constant temperature i.e with fixed ksT = 0.1. N = 12, 8 = 0.8, 9 = 0. 

phase transition at this value of B. The transition becomes progressively blurred at higher 
temperatures and ultimately vanishes at a critical temperature T = T c above which thermal 
fluctuations dominate destroying all pairwise entanglement in the chain. But apart from 
these features also available in the non-dimerized version of Heisenberg chain as observed by 
several authors, there are certain distinctive characteristics available for the case of its dimerized 
version. First of all the values of concurrences in strongly coupled pairs are almost maximal 
(> 0.9) while it is zero for weakly coupled pairs even for 5 = 0.2 (see figs. [2]& [3]) in absence 
of any external magnetic field. Moreover particularly at low temperatures, there exist certain 
other discrete values of B , for each of which the pairwise concurrence for the strongly coupled 
pairs abruptly diminishes to a lower value and then remains steady with the increase of B, 
until next such critical value of B arrives; giving rise to a staircase structure (see fig. [2]), which 
is easier to identify when one takes a constant temperature 2-d slice at a low temperature like 
the figure [5] and [5j It may be apparent from these figures (figs. |4] and [5]) that this staircase 
structure is more pronounced with smaller values of dimer strength S. 



Figure 6: Variation of entanglement within weakly coupled pair with magnetic field at constant 
temperature i.e with fixed ksT = 0.1. N = 12, 5 = 0.2 ,6 = 0. 

The same observations for weakly coupled pairs of the same chain (fig. [3]) also revealed such 
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staircase structure (particularly with small 8, as in fig. [6]) but that structure becomes blurred 
with higher values of 8, except a certain surge near the final B c (see fig. [7]). It may be noted that 
at the extreme value of 8 — 1, the weakly coupled pairs get decoupled and hence entanglement 
within those pairs vanishes at that limit. Moreover it is apparent from the figs. [6] & [7] that 
at zero magnetic field weakly coupled pairs are completely separable and they get entangled 
only after a specific value of external magnetic field sets in (depending upon 5), which may be 
considered as a case of magnetic entanglement. 



Figure 7: Variation of entanglement within weakly coupled pair with magnetic field at constant 
temperature i.e with fixed k B T = 0.1. N — 12,8 — 0.8 ,9 = 0 



Figure 8: Variation of concurrence in strongly coupled pairs with the variation of B and 8 at 
a fixed temperature k B T = 0.1. 

Next the variation of concurrence within strongly coupled pairs with that of external mag¬ 
netic field and that of dimer strength (see figj8]) at a constant low value of temperature shows 
that this staircase feature is more predominant for relatively low values of dimer strengths. 

From the 3-D plot (fig. |8]) it is clear that, if 8 increases from 0 even by a small amount (i.e. 
making it non-dimer chain) the staircase structure of the variation of concurrence with external 


6 




magnetic field appears and the stair sizes vary with 5 and it vanishes when S approaches the 
value 1, the limit in which the chain breaks into a collection of identical spin pairs interacting 
with XXX Hamiltonian. 

The staircase structure may be contributed by the gapped nature of the energy spectra for 
dimer chains, but further investigation is required to identify the actual mechanism. 

To investigate the effect of inclination of external magnetic held, the dimerized XX chain is 
chosen. The variation of the inclination of magnetic held has also been studied and for different 
angles the variations of the pairwise entanglement with held strength are plotted in hgs. BH3 
& 11 It has been found that for XX dimer chain even the strongly coupled pairs develop less 


amount of entanglement as expected intuitively, but it is observed that at low temperature 
with the increase of B , the entanglement hrst diminishes to almost zero at a critical value of 
B and then again build up; moreover this build up increases with the relative strength of the 
longitudinal component of external held, i.e. B x or in other word with the angle of tilt. Similar 
secondary build up of entanglement was also reported by some authors earlier for two qubit 
anisotropic Heisenberg chain with the variation of anisotropy (see (3|fl4)). 



Figure 9: The variation of concurrence in strongly coupled pair [5 = 0.2] in XX model with B 
and ksT for tilted B at an angle 9 = 0 with z-axis. 



Figure 10: The variation of concurrence in strongly coupled pair [<5 = 0.2] in XX model with 
B and ksT for tilted B at an angle 6 = 7t/4 with z-axis. 

The studies already mentioned correspond to closed chains which have translational sym¬ 
metry. But as this study has been motivated by the communication aspect of the dimers as 
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Figure 11: The variation of concurrence in strongly coupled pair [h = 0.2] in XX model with 
B and ksT for tilted B at an angle 6 = n/2 with z-axis i.e. along x-axis. 


recently proposed by some authors (refs), it is natural to extend these studies for open XXX 
dimer chains. Due to absence of translational symmetry each pair has some unique feature in 
phase transition. We have chosen one strongly coupled one at the edge (fig. 12) and another 
one near the middle of the chain (fig. 13). 



Figure 12: Variation of concurrence in the first strongly coupled pair of an open XXX chain 
with external magnetic field B for 5 = 0.2, N = 12 


ft is clearly seen that the staircase structure is retained for the open chains too and in 
absence of external held strongly coupled ones are almost maximally entangled. The pair at 
the edge almost replicates (fig. 12) the features corresponding to strongly coupled pairs in closed 
chains; but for the pair near the middle shows (fig. 13) a build up at some intermediate values 
of B, which is not found for strongly coupled pairs in closed XXX chains. 
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Figure 13: Variation of concurrence in the fifth pair (strongly coupled) of an open XXX chain 
with external magnetic field B for 5 = 0.2, N = 12 
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